1 



Is anisotropic flow really acoustic? 
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The flow harmonics for charged hadrons (v n ) and their ratios (v n /v2) n >3, are studied for a broad 
range of transverse momenta (pr) and centrality (cent) in Pb+Pb collisions at ^/snn = 2.76 TeV. 
They indicate characteristic scaling patterns for viscous damping consistent with the dispersion 
relation for sound propagation in the plasma produced in the collisions. These scaling properties 
are not only a unique signature for anisotropic expansion modulated by viscosity, they provide 
essential constraints for the relaxation time and a distinction between two of the leading models for 
initial eccentricity. These constraints could be important for a more precise determination of the 
specific shear viscosity rj/s. 

6 PACS numbers: 25.75.-q, 25.75.Dw, 25.75.Ld 

7 Azimuthal anisotropy measurements are a key ingre- 
s dient in ongoing efforts to pin down the precise value 

9 of the transport coefficients of the plasma produced in 

10 heavy ion collisions at both the Relativistic Heavy Ion 

11 Collider (RHIC) and the Large Hadron Collider (LHC). 

12 The Fourier coefficients v n are routinely used to quantify 

13 such measurements as a function of collision centrality 

14 (cent) and particle transverse momentum px] 



35 to the magnitude of both e n and the transport coefficient 

36 rj/s (i.e. the specific shear viscosity or ratio of shear vis- 

37 cosity 77 to entropy density s) of the expanding hot matter 

38 [il [Til. Il3 - fl8| . Thus, v n (pT, cent) measurements provide 

39 a crucial bridge to the extraction of rj/s from data 

40 Initial estimates of rj/s from v n measurements 

41 

0, EI US-HI have all indicated a small value (ri/s 
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where <j> is the azimuthal angle of an emitted particle, and 
yj n are the azimuths of the estimated participant event 
planes [JQ; 



v n (pi 



(cosn(0 - ip n )) 



is where the brackets denote averaging over particles and 

19 events. The distribution of the azimuthal angle difference 

20 (A0 = 4> a — 4>b) between particle pairs with transverse 

21 momenta p?p and p\ (respectively) is also commonly used 

22 to quantify the anisotropy 



dN PB.ns 



dAcj) 



l + 5^2t>„,„(i4,i£)cos(nA0) j , (2) 



where the latter factorization has been demonstrated to 
hold well for pt < 3 GcV/c for particle pairs with a 
sizable pseudorapidity gap Arj p Q . 

The coefficients v n (p T ,cent) (for p T < 3 — 4 GeV/c) 
have been attributed to an eccentricity-driven hydrody- 
namic expansion of the plasma produced in the colli- 
sion zone That is, a finite eccentricity moment 
e n drives uneven pressure gradients in- and out of the 
event plane ip n , and the resulting expansion leads to the 
anisotropic flow of particles about this plane. In this 
model framework, the values of v n (pT, cent) are sensitive 



42 times the lower conjectured bound of 1/47T [25|). Recent 

43 3+1D hydrodynamic calculations, which have been quite 

44 successful at reproducing v n (pT, cent) measurements [26l — 

45 [28} , have also indicated a similarly small value of rj/s < 

46 2/47T. However, the precision of all of these extractions 

47 has been hampered by significant theoretical uncertainty, 

48 especially those arising from poor constraints for the ini- 

49 tial eccentricity and the relaxation time. One approach 

50 to the resolution of this issue is to target these uncertain- 

51 ties for systematic study, with the aim of establishing re- 

52 liable upper and lower bounds for 77/s [l2|, H^. An alter- 

53 native approach, adopted in this work, is to ask whether 

54 better constraints for these theoretical bottlenecks can 

56 be developed to aid precision extractions of 77/s? 

57 Given the acoustic nature of anisotropic flow (i.e. it is 

58 driven by pressure gradients), a transparent way to eval- 

59 uate the strength of the dissipative effects which reduce 

60 the magnitude of v n (pT-, cent), is to consider the attenua- 

61 tion of sound waves in the plasma. In the presence of vis- 

62 cosity, sound intensity is exponentially damped e( -r / rs '' 

63 relative to the sound attenuation length T s . This can 

64 be expressed as a perturbation to the energy-momentum 

65 tensor Ta V (3l| : 



6T u „(n,t) = cxp (-/3n 2 )ST u ,(0), f3 



H^J-i (3) 



66 which incorporates the dispersion relation for sound 

67 propagation, as well as the spectrum of initial (t = 0) 

68 perturbations associated with the eccentricity moments. 

69 The latter reflects the collision geometry and its associ- 

70 atcd density driven fluctuations. Here, the viscous coeffi- 

71 cient P oc rj/s, t cx R is the expansion time, T is the tem- 

72 perature, k = n/R is the wave number (i.e. 27ri? = nX 
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FIG. 1. (a)-(d) v„/e„ vs. n for charged hadrons for several pr selections in 20-30% central Pb+Pb collisions at ^/snn = 2.76 

TeV; (e) f3 vs. pr for the same centrality selection. The v n data are taken from Refs. 0, the dashed and dotted curves 
represent fits (see text). 



73 for n > 1) and R is the transverse size of the collision 

74 zone. 

75 The viscous corrections to v n implied in Eq. [3J do not 

76 indicate an explicit pT-dcpendcncc. However, a finite 

77 viscosity in the plasma results in an asymmetry in the 

78 energy-momentum tensor which manifests as a correction 

79 to the local particle distribution (/) at freeze-out (23 



/ = fo + Sf(pr), Pt = -jt; 



(4) 



where /o is the equilibrium distribution and <5/(pr) is 
its first order correction. The latter leads to the px~ 
dependent viscous coefficient ft (pr) oc ft/p T , where the 
magnitude of a is related to the relaxation time tr(j>t)- 
Equations [3] and [4] suggest that for a given central- 
ity, the viscous corrections to the flow harmonics v„ (pt), 
grow exponentially as n 2 ; 

Vnipr) 



oc exp I — p n 



(-ft'r, 



(5) 



87 and the ratios (v n (pT) / V2(pt))u>3 can be expressed as; 
v n (pr) e n _____ ( oV __ 2 ^ 



■ exp 



V 2 (Pt) £2 

88 indicating that they only depend on the eccentricity ra- 

89 tios and the relative viscous correction factors. Note as 

90 well that Eq. [5] shows that the higher order harmonics 
si Vn,n>3, can all be expressed in terms of the lower order 

92 harmonic v 2 , as has been observed recently [1, [32| . For a 

93 given harmonic, Eq. [S]can be linearized to give 



In 



Vn{PT) 



R 



(7) 



which indicates a characteristic system size dependence 
(1/R) of the viscous corrections. 

If validated, the acoustic dissipativc patterns summa- 
rized in Eqs. [5j[6] and [71 indicate that estimates for a, ft 
and £n/£2 can be extracted directly from the data. Here, 
wc perform validation tests for these dissipativc patterns 
with an eye toward more stringent constraints for tr, rj/s 
and the distinction between different eccentricity models. 

The data employed in our analysis are taken from mea- 
surements by the ATLAS collaboration for Pb+Pb col- 
lisions at y/SNN = 2.76 TeV [g], [3(J and p+Pb collisions 
at v/sjvjv = 5.02 TeV [33|. These measurements exploit 
the event plane analysis method (c.f. Eq. [T]), as well as 
the two-particle A</> correlation technique (c.f. Eq. [2]) to 
obtain robust values of v n (pT, cent) for a sizable Ar) p gap 
between particles and the event plane, or particle pairs. 
We divide these values by e„(cent) and plot them as a 
function of n, to make an initial test for viscous damp- 
ing compatible with sound propagation in the plasma 
produced in these collisions. Monte Carlo Glauber (MC- 
Glaubcr) simulations were used to compute the number 
of participants N part (cent) and e„(cent) from the two- 
dimensional profile of the density of sources in the trans- 
verse plane p s (r±). The weight u>(r±) = rj_™ _34_ was 
used to compute e„(cent). 

The open circles in Figs. [1] (a)-(d) show representa- 
tive examples of v n /e n vs. n for several px cuts, for the 
20-30% centrality selection. The dashed curves which 
indicate fits to the data with Eq. [5j confirm the ex- 
pected exponential growth of the viscous corrections to 
V n , as n 2 . The py-dependent viscous coefficients ft (Pt) 
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FIG. 2. v n /v2 vs. pt for several centrality selections for Pb+Pb collisions at ^/sjvjv = 2.76 TeV. The open symbols show the 
values obtained from data; the filled symbols show the results of fits to these ratios with Eq. [6] (see text). 



obtained from these fits, are summarized in Fig. [TJ (c); 
they show the expected dependence attributable to 
5f{pr)- The dotted curve in the figure is a fit which gives 
the values a ~ 0.58 and (3 ~ 0.126. Similar results were 
obtained for a broad range of centrality selections. 

Additional constraints can be obtained from the ratios 
of the flow harmonics (v n (j)T) / V2(pr))n>3 (cf. Eq. as 
well as the dependence of v n (pT)/£ n on the transverse 
size of the collision zone (cf. Eq. [7J . The open symbols 
in Fig. [5] show the values of (v n (pT) / v 2 (pt)) for n = 
3, 4 and 5, for each of the centrality selections indicated. 
A simultaneous fit to these ratios was performed with 
Eq. [5] to extract j3 and e n /E2 at each centrality. Small 
variations about the previously extracted value of a ~ 
0.58 were used to aid the convergence of these fits. The 
filled symbols in Fig. [2] show the excellent fits achieved; 
they confirm the characteristic dependence of the relative 
viscous correction factors expressed in Eq. [5] They also 
confirm that the relationship between i>2 and the higher 
order harmonics stems solely from "acoustic scaling" of 
the viscous corrections to anisotropic flow. The extracted 
values for £ rl /£ 2 , a and f3 are summarized and discussed 
below. 

«9 Figures [3Ja) and (b) gives a more transparent view of 

150 the influence of system size on the viscous corrections. 

151 Fig. |31[a) shows that u 2 ,3 increases for 140 < N part < 340 

152 as would be expected from an increase in £2,3 over the 

153 same N part range. For N part < 140 however, the de- 

154 creasing trend of i>2,3 contrasts with the increasing trends 

155 for £2,3, suggesting that the viscous effects due to much 

156 smaller system sizes, serve to suppress t>2,3- This is con- 

157 firmed by the dashed curves in Fig. [3Jb) which validate 



i the expected linear dependence of ln(v n /e n ) on 1/R (cf. 
1 [7]) for the data shown in Fig. Eta). It is noteworthy that 
1 the p+Pb data point is compatible with this trend. A 
similar dependence was observed for other pt selections. 
: The slopes of these curves serve as an important addi- 
i tional constraint for (3. 

Figures^c) - (e) show a comparison between the £„/£2 
. ratios extracted from the fits shown in Fig. [5] (open sym- 
i bols), and those obtained from model calculations (filled 
■ symbols). For the 5-50% centrality range, the compari- 
i son shows good agreement between the extracted ratios 
1 and those obtained from MC-Glauber calculations with 
1 weight w(r_i_) = r± l 34]. A similarly good agreement 



with the ratios obtained from a Monte Carlo implemen- 
tation 1351 of the factorized Kharzeev-Lcvin-Nardi (KLN) 
model [36j, |37| is not observed. For the 0-5% most cen- 



tral collisions, the extracted values of e n /s2 are signif- 
icantly larger than the values obtained from either ec- 
centricity model. This difference is compatible with an 
additional source of fluctuations in very central events, 
that is not appropriately accounted for in the initial ec- 
centricity models considered. Such fluctuations could be 
an important impediment for precision extraction of 77/s 
from hydrodynamic calculations which emphasize ultra- 
central collisions pjj . 

The fits shown in Fig. [5] also give values for a and 
j3, which are summarized in Figs. [3jf) and (g); they are 
essentially independent of centrality. This suggests that, 
within errors, the full data set for v n (j)T, cent) can be un- 
derstood in terms of the eccentricity moments coupled to 
a single (average) value for a and j3 (respectively). This 
observation is compatible with recent viscous hydrody- 
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FIG. 3. (a) 112,3 vs. Np art for pr = 1 — 2 GeV/c: (b) ln(v n /s n ) vs. 1/.R for the data shown in (a): (c - e) centrality dependence 
of the £„/£2 ratios extracted from fits to (v n (pT) /v2(pt))u>3 with Eq.[6] e n /£2 ratios for the MC-Glauber [ULI^] an d MC-KLN 
[35l ] models are also shown: (f) extracted values of /3 vs. centrality: (g) extracted values of a vs. centrality (see text). 



191 namical calculations which have been successful in repro- 

192 ducing v n (pr, cent) measurements with a single 8f{ftr) 



27| . Therefore, these 



193 ansatz and a single value of 77/s [2f 

194 values of a and (3 should provide an important set of con- 

195 straints for detailed model calculations. 



222 for a more precise determination of the specific shear vis- 

223 cosity 77/s. 
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To demonstrate their utility, we have used one of the 
(3 values to make a rough estimate of Airrj/s ~ 1.5 for 
central collisions for (T) ~ 300 MeV and t ~ 12 fm 
[3^. This estimate, which is approximately 25% larger 
than a similar estimates for central Au+Au collisions at 
RHIC (-y/sjviv = 0.20 TeV) [32} is in reasonable qualita- 



tive agreement with recent estimates [2f 
Detailed model calculations with improved constraints 
for e n and 8f(pr) { t r(pt)) are required to make more 
precise estimates of r\j s. However, our analysis clearly 
demonstrates the value of the relative magnitudes of v n 
for developing such constraints. 

In summary, we have presented a detailed phenomeno- 
logical study of viscous damping of the flow harmonics 
v n and their ratios (v n /(v2))n>3, for Pb+Pb collisions at 
y/sNN = 2.76 TeV. Within a parametrized viscous hydro- 
dynamical framework, this damping can be understood 
to be a consequence of the acoustic nature of anisotropic 
flow. That is, the observed viscous damping reflects the 
detailed scaling properties inferred from the dispersion 
relation for sound propagation in the plasma produced 
in these collisions. These patterns give a unique signa- 
ture for anisotropic expansion modulated by viscosity, 
and provide straightforward constraints for the relaxation 
time and a distinction between two of the leading models 
for initial eccentricity. Such constraints could be crucial 
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